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Abstract 

The electron-positron pair production due to the dynamical Schwinger process in a slowly 
oscillating strong electric field is enhanced by the superposition of a rapidly oscillating weaker 
electric field. A systematic account of the enhancement by the resulting bifrequent field is provided 
for the residual phase space distribution. The enhancement is explained by a severe reduction of 
the suppression in both the tunneling and multiphoton regimes. 
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I. INTRODUCTION 


The developement of new radiation sonrces offers opportnnities for the investigation 
of fnndamental physical processes which hitherto conld not be accessed experimentally. 
In the realm of qnantnm electrodynamics the Schwinger effect, i.e. the decay of a pnre 
electromagnetic held into electron-positron pairs, is among the challenges which have escaped 
an experimental verification nntil now. While pair prodnction via pertnrbative or even 
non-linear effects in different forms, mostly inclnding nnll fields, is well established, the 
genuinely non-perturbative Schwinger effect [1] in a static homogenous electric held requires 
held strengths far beyond current laboratory capabilities. Much hope has been put on optical 
laser facilities where, in fact, high held strengths can be achieved in the focal spots near to 
the refraction limit. However, only future developements as, e.g., the pillar-3 part of ELI [2], 
could bring the dynamical Schwinger effect (which refers to an alternating electric held) into 
realistic reach. Here, synchronized counterpropagating laser beams with suitable polarisation 
can produce regions of an alternating held with a dominating electric component in the 
vicinity of the antinodes. 

In the search for conhgurations which could enable a verihcation of the Schwinger effect 
as a tunneling process the idea has been put forward to enhance the pair creation rate by the 
assistance of a multi-photon process [3, 4]. This set-up is denoted as the assisted dynamical 
Schwinger process. There, one can think of the combination of optical high-intensity laser 
and XFEL beams, the latter ones representing a fast weak held. In fact, at LCLS such 
instrumental prepositions are already at our disposal [5], and the HIBEF collaboration [6] at 
the European XFEL in Hamburg plans an analogous installation, albeit with different key 
parameters w.r.t. repetition rates, energies and intensities. Further experimental proposals 
can be found in [7-9]. Whether such conhgurations will enable to investigate the Schwinger 
effect, or a variant thereof, needs to be elucidated. 

Given such a motivation, we are going to consider here a model for pair production by 
a bifrequent, spatially homogenous electric held which acts for a hnite time interval. The 
pair creation process is a non-equilibrium process to be described by a quantum kinetic 
equation with a strong non-Markovian feature. Such a framework has been employed in 
various previous analyses, e.g., for the superposition of Sauter pulses [10, 11] or periodic 
helds and Sauter pulses [12] or for two (or more) periodic helds [13, 14], and for studying 
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the temporal structure of particle creation [15-17]. Alternative frameworks make use of 
WKB-type approximations [18], worldline instantons [19, 20] or lightfront methods [21, 22], 
Our goal is to investigate the dynamical Schwinger effect in bifrequent helds systematically 
over a large region in the parameter space spanned by held strengths and frequencies. 

II. RESIDUAL PHASE SPACE DISTRIBUTION: ANALYTICAL APPROXIMA¬ 
TIONS 

The residual phase space distribution / of pairs in a bifrequent electric E{t) held (fre¬ 
quencies z/i, z /2 = Nvi) which acted for a hnite time ff.t. with constant amplitudes Ei and E 2 
is given by (see Appendix) 

^ 00 

e>imin 

with momentum components perpendicular {pj_ = p cos <p) and parallel (py = p sin ip) to E 
(p = —I ■ ■ ■ |). The function S{x) = sin^ x/x^ has a main peak at a; = 0, S'(O) = 1, of width 
27r and an inhnite number of side peaks declining as ~ x~‘^. This implies that the phase 
space distribution (1) for hxed p displays a series of main peaks at p = pe{p), see Fig. 1. 
(Accounting for the p degree of freedom gives rise to ridges over the p±-p\\ plane, sometimes 
called shell structure.) The peak positions pe{p) are determined by the resonance condition 
(cf. |23|) 

2Sl(p,(ip), ip) - (2) 

27V 1 

^(p±,P||) = ^yda;|^l + ^+(^^-7f^cosa;-72"^cosiVa;) ^ , (3) 

0 

where 71,2 = Keldysh parameters with Ec = —e the electron charge and m 

its mass. Since Q{p£{p), p) > > m with the effective mass m* = fl(p£ = 0) (cf. [24, 25] 

for the ehective mass concept in the single held case), the values of ^ must exceed a certain 
threshold, £ > where [a;] denotes the smallest integer greater than or equal to 

X. 

The number of side peaks belonging to pe until the next main peak at p^+i is about 
k = (p£+i —pe) ^'{pe{p), p) tf.t./TT 3> 1. This means that the spectrum consists, for sufficiently 
large ff.t,, of well separated main peaks centered at pe(p), with some micro-structures emerging 
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Figure 1: Residual phase space distributions at py = 0 (solid curves, from integrating the quantum 
kinetic equations, e.g. equations (1-7) in [23]) for field “1” (left panel, Ei = 0.2Ec, i^i = 0.02m, 
i.e. 7 i = 0.1), fields “l-|-2” (middle panel) and field “2” (right panel, E 2 = O.ObEc, 1^2 = 0.5m, i.e. 
72 = 10). The dashed curves depict the spectral envelope £{p±) with S given by (5). The 

flat-top duration is determined as in [23] by = 27r • 50 and the C°° ramping (switching-on) 

and deramping (switching-off) periods are such to cover 5 oscillations each. 



from the superposition of the side peaks encoded in S, which are not resolved on the scale 
displayed in Fig. 1 . Due to lim S{atf,t.) = ‘2TTtjl5{a), Eq. ( 1 ) asymptotically gets a form 

n.t.-s-oo 

derived in [26] (see also [27]), albeit for a single field and in a low-momentum expansion and 
for the rate, which is time (tf.t.) independent. We consider here hnite values of tf.t., (i) on the 
one hand sufficiently large to accomodate at least a few oscillations of both fields “ 1 ” and “ 2 ” 
within the flat-top (tf.t.) period of the overall shape function of the potential, and (ii) not 
too large to avoid the onset of Rabi oscillations [14] which would modify the tf.t. dependence 
in ( 1 ). 

The main peak altitudes are governed by the squared Fourier-coefficients \F£[p£{(p), ip) 
which contain the spectral envelope function S{p±). As pointed out in the Appendix, for 
P|l = 0 , they can be approximated by 


\Fe{pi,p = 0)1^ ^ -iy^S{pe){l - (-l)'^) , 
S{p^) = ^ 


XQ 


G(p±, 71 , 72 ,^) = jdx ( 1 + ^ - ( 7 r^ sinhx-f 72 “^ sinhiVa;)' 


( 4 ) 

( 5 ) 

( 6 ) 
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for the special case iV = 4n + 1, n = 0,1, 2, • ■ ■. The quantity xq is determined by 

7 ]“^ sinhxo + 7 ^^ sinhiVxo = \Jl + p\/m?‘. (7) 

In the single held case {E 2 —t 0, 72 —)■ 00 ), G reduces to 

Gi{p±,-fi) = G(72 ^ 00 ) = ^7i (^1 + ^^g 1^71^1 + , (8) 

where g{'y) = iT 2 (|, 2, — 7 ^) is a function already introduced in [24], cf. also [26]. We argue 

that a handy approximation is provided by 



and, analogously. 



The former is appropriate for xq < 1 and the latter for \l ^ ^ li < 10. The spectral 

envelopes are displayed in Fig. 1 by dashed curves using ( 6 ). Evidently, the envelopes connect 
the peaks of the fairly rich spectra very well, the details of which are only accessible by 
integrating the full set of quantum kinetic equations. The approximations (9) and (10) are 
suitable for conditions as in the left panel of Fig. 1, while (9) gives some semi-quantitative 
account (within factor of two) for the middle panel. 

To understand qnalitatively the amplihcation by a second weak held, as exemplihed 
in Fig. 1 (fnrther examples in [23]), one can nse the approximation (9) for p± m and 
the hrst-order iterative solntion of (7), Xq = arsinh ^71 — ^ sinh[iV arsinh 7 i]^. Snpposing 
7 i 1 and iV 7 i 1 one hnds for large 72 

G(7i,72,iV)^Gi(7i)--lV, (11) 

72 

i.e. dne to the presence of the second held, encoded in 72 , the modnlus of the exponential, 
4^Gi, becomes diminished by abont 4 ^ interestingly independent of z /2 in the 

given approximation (which is a special case of the expansion of a mnlti-scale implicit 
fnnction). In other words, the spectral envelope of held “1” gets lifted by exp 
dne to the presence of held “2”. In general we stress that, dne to the monotonic behavior 
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of the hyperbolic sine functions in ( 6 ) and the dehning equation (7) for xq, one infers 
a;o(7i,72 —)■ oo, TV) > 0:0(71,72,^) > 0 and analogously Gi(7i) > (7(71,72, TV) > 0. Since 
the negative of G enters the exponential of the envelope function, a dropping of G by the 
second weak held causes the anticipated amplihcation. A quantitative consideration is given 
in the next section. 

III. DISCUSSION OF THE AMPLIFICATION 

A. Single field case 

In the following we consider £ = exp 4()^(7(px, 7i, 72, A) j as the essential part of 
the spectral envelope of the phase space distribution. The benchmark for the following 
discussion is the reference distribution for one periodic held, i.e. 72 —)■ 00 as realized by 
B 2 —)■ 0. The exponential 4^Gi of the spectral envelope function Ti is displayed in Fig. 2. 
The region 71 < 1 is often termed adiabatic (tunneling) region where the residual pair 
density is independent of frequency, while 71 > 1 is the anti-adiabatic (multi-photon) region 
which is weakly dependent on the held strength. Focusing on low-momentum particles, 
Px/m 1, the spectral envelope in the adiabatic region behaves as exp(— 7r(i7c/77i)), while 
in the anti-adiabatic region it behaves like exp(—4^ log 471 -|- 4^) in leading order based 
on (9). The contours (solid curves) in Fig. 2 are based on ( 6 ); they differ marginally on the 
approximate estimates (9) (dotted curves). Figure 2 provides a quantitative description of 
the landscape of the dynamical Schwinger effect, hrst qualitatively discussed in [15]. Note 
the huge variation over the displayed parameter range by recollecting that the residual phase 
density at py = 0 is bounded from above by with Ti = exp 4^(7ij, i.e. in the 

blue region one meets a very strong suppression of the pair density, guaranteeing a long 
lifetime of the vacuum. Only the red region above the 71 = 1 line seems to allow for verifying 
experimentally Schwinger’s tunneling process in one oscillating held. As pointed out in [23] 
(see also [14] and further references therein), the impact of a second held can signihcantly 
reduce the strong suppression due to large values of 4^G, i.e., it results in an amplihcation 
ehect, as mentioned above. 
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z/i/m 


Figure 2: Contour plot of the exponential 4^Gi = 4:^G{p± <C m-, 71,72 —)• 00 ) of the single-held 
spectral envelope function over the held-frequency {Ei/Ec vs. ui/m) plane (solid curves: using ( 6 ) 
for G, dotted curves: the approximations (9,10)). Note that the spectral envelope is given by 
£i with £i = exp 4^Gi^. Lines (dashed) of constant Keldysh parameters are displayed 
too. Parameters from table 1 in [27] are depicted by stars, while the bullet is for ELI-NP [28]; the 
ellipse is a former perspective of ELI used, e.g., in [3, 4]. The triangle is for MaRIE at LANL [29]. 
Point 1 and the square 2 depict parameters considered in [23], where the validity of Eqs. (l)-( 6 ) 
can be checked by numerical integration of the quantum kinetic equations, as is also exemplihed in 
Eig. 1. Light grey dashed contour curves are for 2, • • • ,9 between the solid decade contour curves. 
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B. Two field case: Amplification 


The key for an interpretation of the amplification effect by an assisting field (z/2, 
E 2 ) is the approximation (9) for p_i_ m: Since 0:0(71,72,-/V) < 0:0(71,72 00 ), 

the spectral envelope fnnction Si = exp 4^Gij is lifted by the amonnt S 1 + 2 /S 1 = 
exp(4^[-G'(7i,72,iV) + G'i( 7 i)]) exp(7r^ [-0:0(71,72, + a:o(7i, 72 ^ 00)]), where 

Ti +2 = exp 4^G(7 i, 72, iV) j. Dne to the nsnally considered range of valnes even a 
moderate reduction of G (or 0:0) due to the assisting field, e.g. by (9(50%), leads to a huge 
reduction of the suppression in the subcritical region Ei/E^. ^ 1 and vi/m 1 steered either 
by '^'yi = Ec/El in the adiabatic region or by log 71 in the anti-adiabatic region. A few 
examples are exhibited in Fig. 3, where we show the exponential 4^(9 (p_l m,7i,72,iV), 

entering the spectral envelope |z/ft£^ Ti+2, for a few selected fields i?i, ui over the E 2 -i '2 
plane. As reference one has to take the spectral envelope function from Fig. 2. The value 
of the envelope exponential is trivially reduced by 50% in the points E 2 = Ei and z/2 = vi 
(marked by symbols in Fig. 2), i.e. simply doubling the field strength. This is the lowest order 
of the multi-field configurations considered in [30]. Further amplification occurs for enlarging 
z/2, as signalled by a reduction of 4^(9. That is, the amplification beyond the field doubling 
happens in the region right to the contour curve (heavy grey dashed) going through the 
reference point (bullet) E 2 = Ei, z/2 = z/i. The single field result is recovered by inspection 
of E 2 —)■ 0, i.e. going down from the reference point at z/i, Ei. Note that our amplification 
estimate by formula (4) requires E 2 < Ei and z/2 = Nui with N = 4n-|-l, n = 0,1,2---, 
that is our result holds true actually only at discrete values of z/2. 

The net outcome points to the importance of the frequency z/2, while variations of E 2 
are subleading. To achieve amplifications which overcome the strong suppression in the 
tunneling regime one has to employ frequencies z/2 near or above critical range z/ > m, 
where essentially the multi-photon process sets the scale, i.e. one has to compare the phase 
space population with the one for the z/2, E 2 field alone. This is accomplished by inspecting 
Fig. 2 for the exponential of the spectral envelope function and by correcting the pre¬ 
exponential factor in (4) by z/i —)■ z/2. As shown in [23] for selected examples (see also 
Fig. 1), the action of the field (z/2, E 2 ) lifts the phase space distribution emerging from (z/i, 
El) alone above the one of (z/2, E 2 ) alone. A decisive difference is in the density of states 
(peaks) in the p±{p\\) distributions; The number of peaks within pj_ = 0 - ■ - m is given by 



o 


o 


CO 

O 


o 



^ 3/^3 


^ 3/^3 


o; 

'a 

OP 

t:3 

c/2 

o 


o 

s 

o 

’bb 

OP 

?H 


-l-i 

s 

a 

OP 

-l-i 


tq 


CO 

2 

2 

ca 

'bb 

*52 


CM 

c^ 


V 

H 

a, 

o 

B\ ^ 

2 

*-+^ 

a;> 

o 

a 

>< 

a; 

OP 


o 


?H 

o 

o 

U 


a:> 

?H 

.£P 

2 






V 

H 

SI i 

r-^ 

2 

-l-i 

03 


02 

cd 




02 

bb 

bjO 


CM 


+ 


TU 

o3 

V 

r-'^ 

tq 

fH 

2 

2 


bJD 

o 


o 

bO 


T3 

02 

02 

Sh 

CO 

O 

02 

02 

cC 

CO 

?-H 

02 

T3 

02 

CO 

t3 

02 

bO 

02 

02 

H 


CO 

Sh 


CO 

02 

> 

O 

02 

ri:^ 

CO 

o3 

XJ 

02 

U 

bJD 

02 

02 


bC 


o 

'bb 

02 

xi 

02 

O 





o 

02 

02 

> 

r-^ 

■+^ 

?H 

o 

02 

02 

O' 

_o 

o 


02 


02 

-S 

, Sh 

1 

CO 

02 

T3 

CO 


02 

cC 


"02 


02 

02 


o3 

02 

02 

02 

CO 

T3 

O 




bC 


3 


-S 

i-t 

02 

CO 

02 

CO 


> 


02 


O 

'B 

O 


O 

"bO 

-S 

02 




02 



O 


Of^ 


'+J 

-4-i 

02 

o 


02 

t:! 

:2 



"02 

CO 


t:J 




02 

«-M 

Q 





c^T 



O 


02 

O 

a 

?-H 

.2 

02 

CO 

02 

'S 

02 


TU 


-4-^ 


og 

02 

;-i 

03 


CO 

02 

> 


'ft 

!-h 

g 

<jj 

CO 

02 

03 

‘rS 

02 

?-H 


> 


<i 


.r:; 

02 


02 

Q 


o 

02 

T3 

02 

CO 

cC 

T3 


.r:; 

02 


9 




















































Figure 4: Contour plot of the number of main peaks in the inverval p_L = 0 • • • m at py =0 for 
the action of one field characterized by v and E. Light grey dashed contour curves are for 2, • • • ,9 
between the solid decade contour curves. Dashed lines depict constant Keldysh parameters. 

^ = m,ip = 0, 7 i, 72 , N) — n{p = 0,ip = 0, 71 , 72 , iV)] , which can be much larger for the 

field (z/i, El) alone than the field (z/2, E 2 ) alone, depending on the locus in the v-E plane, 
see Fig. 4. 


IV. DISCUSSION 

The goal of our paper is to provide a qualitative understanding for the amplification of the 
pair creation rate in a periodic electric field under the resonant assistance of a second faster 
and weaker electric field. In the quantum kinetic approach, such an investigation would 
require scans of the four-dimensional parameter space vi^ for a dense set of points in 
the two-dimensional pi.-p\\ phase space. In addition, the impact of the temporal structure 
(details of switching on/off and duration of the flat top time span) should be considered. 
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To avoid such a cumbersome wealth of information we refrain from considering the harsh 
landscape of the phase space distribution by analysing only its spectral envelope instead. As 
a consequence, we do not need to consider the details of the peak (ridge) positions, which 
are modihed by the assistance held. We also rely on the previous observation [23] concerning 
the dependence of the residual phase space distribution on the time structure. This in turn 
is based on a harmonic analysis in the low-density approximation [23]. 

Under these conditions, the amplihcation, observed in several previous investigations which 
are restricted to small patches in the parameter space, can be explained as a reduction of 
the huge suppression provided by exp(—7r(Uc/Ui)a;o) in the adiabatic or by exp(—4(m/z/)a;o) 
in the anti-adiabatic regimes: by the assistance held, the quantity xo (the zero of a simple 
function) becomes smaller. Applying these hndings to typical parameters which represent 
the optical laser-XFEL combination according to [27] we hnd a promising perspective for 
laser intensities only signihcantly above the present ELI-NP plan [28]. More promising is 
the laser-7 beam combination, as already pointed out in [3, 4] for ultra-intense laser beams, 
but we see also a good discovery potential for ELI-NP [28] and even multi-PW lasers. The 
potentially dangerous 7 conversion processes might be rejected by the phase space distribution 
of the residual pairs: The amplihed tunneling production displays the distinct peak (ridge) 
structure dictated by the frequency z/i. 

To arrive at such a simple picture we made two restrictions: {i) N = 4n -|-1, n = 0,1, 2, • ■ ■ 
was selected since the main pole in the complex time plane is then shifted in purely imaginary 
direction towards the real axis under the action of the assistance held (see Appendix), thus 
allowing for an easy access to the spectral envelope function and (ii) p\\ = 0 was chosen as it 
was shown in [23] to provide a proper orientation. The avenues for next generalizations are 
obvious: extensions to arbitrary integer N (and further also real N which go beyond the 
resonant type patterns) as well as further poles and estimates of the phase space-integrated 
residual distribution to arrive at a density. Still, the model is fairly simple and needs further 
considerations of important effects, such as spatial inhomogenities which are, for example, 
elaborated in [20, 31, 32]. 
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V. SUMMARY 


In summary we provide a comprehensive tool to access the amplihcation of the dynamical 
Schwinger process by a resonantly assisting second held, both acting for a hnite duration. A 
promising perspective is seen in the combination of ultra-intense optical laser beams with 7 
beams, while the optical laser-XFEL combination also allows for huge amplihcation effects, 
however, probably not sufficient to enable an experimental verihcation. 
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Appendix: Derivation of Eqs. (l)-(6) 


In the low-density approximation of the quantum kinetic equation [33, 34], the electron- 
positron pair density / is given by f{p,t) = ^ \I{p,t)f with 


Hp,t) 


t 


0 


eE(t')e± 

uj{p,t')^ 




( 12 ) 


where we introduced e_L = = ym‘^+p‘]_ + {p\\ —eA{t)Y^ ^{pA) = 

t 

fdt'u{p,t'). The vector potential A(t) and hence the electric field E{t) = —A(t) are 
0 

assumed to be periodic in time, e.g. A = (Ei/z/i) cos(z/it) -|- {E 2 /{Nh'i)) cos{Nh'it) over a 
finite time span tu., much longer than the switching-on/off intervals. Accordingly one can 
split 0 in a linearly growing and a periodic part, 0(p,f) = Q{p)t + P{p,t), to arrive at 
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Figure 5: The integration path Fi from —Tl2 to +Tl2 along the real axis in (15) is substituted 
by the sequence r2 + r3 + r4 + r5 in the upper half of the complex t plane. Dots depict the zeros 
of where branch cuts (dashed lines) emerge. 


i 

I{p,t) = Jdt'F(p, with F{p,t) = The function F is time periodic 

and is decomposed into Fonrier modes, F{p,t) = F£(p)e“*^'^^ so that / turns into 


HP,t) = ^'iFe{p)- 


iu — 2fl{p) 


(13) 


The quantity I becomes large whenever 2Q{p) — iu is small, yielding (2) and (3). Employing 
polar coordinates one defines the ridge radins Pi{(p) by 2n(^pi{ip), (p) — iu = 0. Since fl(p, p) 
as a function of p is strictly monotonous, this resonance equation has solutions only for 


^ f^min 


20(0,y) 


Expanding every summand in I around its corresponding peak position 


at Pi to first order, and keeping Fi fixed, we arrive at 


Tf r \ nGxp [2ffi'(p£((^),(y9)(p-p^((p))t] -1 

/(P.o= E ^ 


(14) 


£=i„ 


terms with i < do not contribute and are thus dropped. The prime denotes the derivative 
w.r.t. p, i.e. fl'(p, p) = dQ{p,p)/dp. To arrive at / one must take the squared modulus 
of (14). Terms mixing different is in this square will go to zero upon switching off, as one 
can argue from a slowly varying envelope approximation. So only the non-mixing terms 
snrvive for t > tf.t., yielding in particular the sir? x/x^ part of Eq. (1), which contains ridge 
positions and ridge widths as anticipated in [23]. For larger times Ux. and small momenta, 
the Popov formnla [26, 27]) is recovered for the special case E 2 = 0. 

What remains is a formula for the Fourier coefficients F7, which determine the ridge 
heights. This can be achieved by deforming the integration contour in the complex time 
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plane and using the method of steepest descent, similar to [24], For the single held case, the 
Fourier coefficients are given by 

r/2 




eE{t)e± 


-T/2 






(15) 


The function u{ty has four zeros of hrst order in the strip —T/2 < Ret < T/2 in the 
complex time plane at to, to, —to and — tg. We deform the integration contour to the 
sequence F 2 + F 3 + F 4 + F 5 shown in Fig. 5. The contributions from F 2 and F 5 cancel due to 
the periodicity of the integrand. The contribution from F 4 is the negative of the complex 
conjugate of that from F 3 , in symbolic notation = —2i Im To evaluate 

the latter integral approximately, we note that the exponent in the integrand is stationary at 
to- One chooses F 3 such that Im© grows rapidly on F 3 away from to, so that quickly 

goes to zero (method of steepest descent). Then only contributions to the integral close to to 
matter, thus enabling an expansion uj{ty = 2ieo_eE{to){t — to) + (9((t — to)^) yielding (we 
suppress the momentum arguments) 
eE{t)e± 1 1 


u{ty^ 


2it — to 


0 ( 1 ) 


(16) 


The expansion of the exponent uses 0(t) = f*dt'cj(t') = ff^°dt'Lo(t') + jydt'u{t') = ©(to) + 
J^dt'^y2eE{to)e±^yi{t' — to) + C>((t — to)^) with the result 

©(t) = ©(to) - ‘^i^/2eE{to)e± {i{t - to)) " + 0((t - to)5) . (17) 


Inserting (16) and (17) into (15) yields our desired result 


Ee = --ilm 


fdte2*®ho) 


g| V^2eT(t/)7/(*h-io))2 


i{t - to) 


3T 


. (18) 


To generalize to two helds we need to sum over all complex zeros ti of with 0 < Re t* < T/2, 
T = 27r/z/i, and Imt* > 0, leading to 


E^{pi{if),^p) ^ -tz/. (19) 

3 

i 

It is instructive to take (p = 0, that is we consider the phase space distribution for p\\ = 0. 
For the case N = 4n + 1, integer n, a sequence of iV + 1 zeros appears at Imt > 0, 
—T/2 < Ret < T/2 (see Fig. 6 , red squares). In particular, the zeros at ±T/4 + ixo/vi get 
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Imt 


T/IO 





-r/4 


r/4 



Rei 


Figure 6: The distribution of the zeros of in the upper half plane for field “1+2” (red squares) 
for the example of = 25. The two zeros for field “1” are shown as black dots, which get shifted 
down (arrows) upon the impact of field “2”. 

shifted down, since xq is smaller in the two-held case than in the single-held case. The sum 
in (19) is dominated by the contribution with the smallest imaginary part, which is the zero 
to = T/4 + ixo/vi- Keeping this leading term and using furthermore 0(T) = 40(T/4) yields 

^TT Vn f / T)^ ( 1 1 \ ^ 

0(^o) = “r+*— dx\l-\ ——sinhx-| -sinhiVa; , (20) 

4 z/ J y VTi 72 / 

0 

hnally leading to (4) and (6). The amplihcation effect by a second assistance held is thus 
rooted in a shift of the leading-order pole towards the real axis in the complex t plane. The 
condition E 2 < Ei ensures the pattern exhibited in Fig. 6 and iV = 4n + 1 is specihc for the 
shift of the single held zeros (black dots in Fig. 6) parallel to the imaginary t axis. 
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